Collective oscillations driven by correlation in the nonlinear optical regime 
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We present an analytical and numerical study of the coherent exciton polarization including exciton- 
exciton correlation. The time evolution after excitation with ultrashort optical pulses can be divided 
into a slowly varying polarization component and novel ultrafast collective modes. The frequency and 
damping of the collective modes are determined by the high-frequency properties of the retarded two- 
exciton correlation function, which includes Coulomb effects beyond the mean-field approximation. 
The overall time evolution depends on the low-frequency spectral behavior. The collective mode, 
well separated from the slower coherent density evolution, manifests itself in the coherent emission 
of a resonantly excited excitonic system, as demonstrated numerically. 
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Observations of coherent optical phenomena require 
a sufficiently long dephasing time. Thus, coherent pro- 
cesses such as the Rabi oscillations W and form-invariant 
pulse propagation ||] (also known as self-induced trans- 
parency) are well studied both theoretically and exper- 
imentally in atomic vapors [0. In solid-state physics, 
the nonlinear optical properties of semiconductors are 
expected to exhibit analogous behavior, especially under 
resonant excitation of excitons plpl, where the exciton 
transition approximates an ideal two-level system. How- 
ever, there is a distinct difference in semiconductors. At 
moderate excitation densities, the mobile, large-radius 
Wannier excitons suffer strong Coulomb interaction on 
close approach and, at high densities, the collection of 
excitons dissociate into a neutral plasma. 

A mean-field theory of the exciton-exciton interaction 
leads to a purely coherent repulsive interaction, which 
does not yield a mechanism for dephasing due to the 
fluctuating effective local field. This level of approxima- 
tion allows for nonlinear electric field-driven density os- 
cillations, which are termed as Rabi oscillations [§-g|. A 
further low-density approximate relation for the density 
and the polarization leads to a Gross-Pitaevskii equa- 
tion pJlOj] for the exciton dynamics. This yields form- 
invariant soliton solutions fln|-[l3|l . 

Keldysh and Kozlov [jTij included correlation in a low- 
density ensemble of excitons. Here we present an ex- 
tension of the Gross-Pitaevskii equation to include the 
correlation effects and explore some fundamental conse- 
quences in the ultrafast nonlinear optical processes. The 
correlation effects are based on our previous theory of 
exciton correlation in the weak nonlinear optical regime 
and are contained entirely in a retarded memory function 
|Eq,16[. The reactance part of the correlation function 
provides a renormalization of the coherent interaction 
and the admittance part provides a dephasing mechanism 
from exciton-exciton interaction. The latter contribu- 
tion may be regarded as a nonlocal temporal excitation- 



induced dephasing |l7| . A new finding is the ultrafast 
collective density oscillations riding on the slow coherent 
density, both being excited by a short laser pulse. That 
this oscillation survives the excitation pulse shows that it 
is not related to the Rabi oscillations. Because of the in- 
clusion of the correlation effects, which, related by a sum 
rule to the exciton local-field interaction, must be just as 
strong, the collective oscillations and the slow coherent 
dynamics together are expected to give a more accurate 
picture of interacting excitons in the density regime be- 
low dissociation. 

The coherent dynamics of the exciton Jl5| can be de- 
scribed by a nonlinear equation of motion for the polar- 
ization density P of the form: 
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The equation is easily generalized to include the spin, 
higher exciton levels, and the continuum and to account 
for spatial variation. It is then exact to the third order in 
the exciting field [p~5| , p~8|JlS| l - The whole equation is thus 
valid for the low-density regime. Our purpose here is to 
explore the temporal variation of the exciton dynamics 
at low excitation density and, therefore, we retain only 
the resonantly excited exciton transition and set P to be 
the usual polarization density divided by /i|0i s (O)| 2 , /i 
being the dipole moment and <fii s (r) the relative motion 
wave function. For the terms in order on the right side 
of the equation, S denotes the detuning of the central 
laser frequency from the resonance, (3 > is a measure 
of the mean-field exciton-exciton interaction, and n c the 
maximum density divided by |0i s (O)| 2 for the validity 
limit in the Pauli blocking term. The Rabi frequency 
fi is given by the dipole energy of the laser field am- 
plitude. Exciton-exciton correlation is contained in the 
retarded memory function F(t). Similar effective equa- 



tions without the correlation term F were already used in 
the literature H^O 21 1. A simple way to account for the 
correlation effect by connecting the polarization equation 
to an antiparallel-spin biexciton density |22j] is equivalent 
to setting F to be the single-frequency biexciton propa- 
gator. 

The retarded memory function is given by p5| , 



F(t) = V\4> ls (0)\ 2 (0\De- tHT rf\0) 9(r) , 
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in terms of the Hamiltonian of the system H, volume V, 
and the interaction operator between two excitons, 
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+ v(r 2 - r' 2 ) - v( n - r' 2 ) - v(r 2 - ri)]B(/ lt r' 2 ), (3) 

with the Coulomb interaction v(r) and the exciton anni- 
hilator 



B(ri,r 2 ) =^i s (ri - r 2 )i>\{ri)^2{ri), 
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where ip\, ip 2 are, respectively, the electron creator in the 
valence band and the annihilator in the conduction band. 
The formulas in configuration space show clearly the in- 
clusion of the electron-hole attraction in each exciton, 
leaving only the interaction between the constituents of 
the two excitons in the interaction operator D. The cor- 
relation effects contained in the memory function are in 
second or higher order of the interaction between two 
excitons. 

Our numerical simulation below shows that intrinsic 
dephasing times due to exciton-exciton scatterings are of 
the order of a few picoseconds in a GaAs quantum well, of 
the same order as that due to electron-phonon ||. Thus, 
the conclusions drawn from the exclusion of the extrin- 
sic dephasing would not be qualitatively changed by its 
inclusion. The equation for the polarization density re- 
sembles the Gross-Pitaevskii equation for a single com- 
ponent Bose-Einstein condensate p,p!o[. However, there 
are two essential differences, (1) the exciton coherence 
is driven externally, and (2) the effective interaction be- 
tween polarization includes the correlation term beyond 
the mean-field theory, which may well be applied to the 
trapped atoms as their density increases. 

After excitation with a short laser pulse at time t = 0, 
the coherent polarization evolves without the driving field 
il. We present first an analytic solution of Eq. (Ill) under 
the assumption of a sufficient slowly varying polarization 
Po(t), which is then used to understand the behavior of 
the numerical solutions without this approximation. Un- 
der the slowly varying assumption, the rapidly decaying 
memory function can be decoupled from the polarization 



F(t-t')P$(t') 2 dt'nP$(t) F(t-r)dT. (5) 



This is similar to a Markov approximation for the polar- 
ization but the scattering rate is still a time-dependent 
memory function, which we group with the mean-field 
term from Eq. (TO) in a new function: 



g(t)= / F{r)dT + i(3. 



(6) 



After the exciting pulse dies down, the detuning 5 term 
can be included in Pq and removed from the equation 
of motion. For the resultant equation, the ansatz Pq = 
n i/2 e i<f> so i ves t ne a \ ow time dynamics: 



t ^ - 1 



n(t) = n 1 + n / gii(t)dt 



<Kt) = <h-2j 9i(t>(t')dt' 



(7) 



for the initial conditions n(0) = no, 0(0) = cf>o set by 
the exciting pulse. The real part of g(t), gii(t), deter- 
mines the dephasing of the coherent polarization and, 
thus, the exciton population. The imaginary part, gi{t) 
determines the phase of the polarization. The mean-field 
effects are included in the total memory function g(t). 

The normalized spectral density of the D-D correlation 
function, denoted by p(oj), is 1/tt times the real part of 
the Fourier transform of F(t)/F(0). The sum rule p3|, 



F(0) / w~ 1 p{w)du = (3 
Jo 
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shows important cancellation in g(t) between the mean- 
field exchange and the correlation effects. The low- 
frequency behavior of the spectral density is governed by 
a power law, p(u>) ~ U) a , discussed below. It determines 
the asymptotic long-time dynamics via: 
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FIG. 1. Equal-spin correlation function (solid-line) for a 
one- dimensional model system. The inset shows the exact 
spectral density p(uj)/u> for small u, indicating a linear de- 
pendence on u. 



where T is the usual Gamma function, leading to a 
crossover from a decaying coherence for a < 1 to a con- 
stant density n(oo) < n for a > 1. In Figure [l], the spec- 
tral density of the equal-spin correlation function is plot- 
ted for a one-dimensional system defined in Jig ] with the 
exciton binding energy u) x — 6.7 meV, the one-electron 
bandwidth of 50 meV and equal electron and hole mass. 
An artificial broadening of 0.4 meV is used to smooth 
the energy splittings due to the finite simulation size of 
N = 320 sites. The inset shows the exact behavior of the 
spectral density for small w, indicating a linear behavior, 
i.e., a = 1. 
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FIG. 2. Coherent exciton density for the one-dimensional 
model of Fig. fil for a pulse intensity FWHM of 200 fs and in- 
creasing pulse area O = J Q(T)dr. The time unit T x — 2ir/u: x . 



Figure H shows the result of the numerical solution of 
Eq. (|l]) for the coherent evolution of the exciton den- 
sity for the model in Figure [fl, using the derived val- 
ues of 13 = 6.6 meV, F(0) = 48 meV and n c = 0.38. 
The transient time evolution under the excitation with a 
200 fs pulse (intensity FWHM) of increasing area shows 
a hint of density oscillations. The pronounced decay of 
the population shortly after the excitation is due to the 
decoherence of the exciton-exciton interaction, following 
by a weak decay asymptotic behavior of the a = 1 case 
discussed above. This regime will also be affected by de- 
phasing of the exciton due to additional degrees of free- 
dom not included in the model. Accurate correlation 
functions for realistic quantum- wells or bulk systems are 
not yet available. We construct a model correlation func- 
tion for two equal-spin excitons in a two-parabolic-band 
model for a three-dimensional semiconductor. In systems 
of higher dimensions than one, we expect, by a phase 
space argument for the scattering of two excitons, that 
the low-frequency dependence of the spectral density is 
super-linear. Combining a low-frequency power law with 
a high-frequency decay approximated as exponential, we 
arrive at a simple model for the spectral density 



T(a- 



iK +1 
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for positive frequencies with a single parameter a{> 0). 
The values for = 52uj x /3, l/n c = 7, and F(0) « 14u/* 
are computed from the two-band model. The frequency 
scale ujp — F(0)//3a is fixed by the sum-rule M). 
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FIG. 3. Coherent polarization for the three-dimensional 
semiconductor model with uj x — 4.2 meV, pulse width 200 fs 
and pulse area 47r for various a of the model correlation 
function. The other parameters are given in the text. The 
dashed-dotted line is the excitation pulse. 

Figure y shows the numerical solution of Eq. (|l|) for 
the model correlation with various values for the param- 
eter a for an excitation with a 200 fs laser pulse and area 
4-7T. For a > 2, the asymptotic behavior of the exciton 
density changes to constant from decaying. For a > 1, 
small oscillations are superimposed on the slow density 
evolution. The period of oscillations is about 200 fs for 
a = 4. These oscillations do not correspond to the fre- 
quency peak io ' c — F(0)//3 of the spectral density, which 
does not vary with increasing a as the peak narrows. 
They are collective modes in density fluctuation excited 
by the strong 47r-pulse. They do not occur in the slowly 
varying approximation dq). It is, however, possible to 
treat the modes as a small, fast-varying perturbation Pi 
to the slow component Pq. Inserting the combination 
P = Pq + Pi in Eq. p) and keeping Pi only to first 
order, we find 



Pi 



\Po\ 



i(3Pi+ / F{t-t')Pi(t')dt' 



(11) 



^9{t)P*i{t). 



The last contribution of Eq. (p]), being P P*/Pq is neg- 
ligible compared with Pi. The resulting linear equation 
can be solved by Fourier transform of Pi (t) for a slowly- 
varying density to yield the collective-mode frequencies 
as zeros of the equation 



w = nE(w), 
E(w) = - F(0) 



p(u)')dui' 
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Figure [| shows a graphical solution of the collective 
mode equation from the crossings between the real part 
of £ and the straight line to/n. The lo = solution cor- 
responds to the Goldstone mode E3] . Models for a range 
of a parameters are shown. Above a critical density, e.g., 
n ss 0.1 (in the density units of Fig. M) for a — 2, addi- 
tional solutions with two high frequencies exist, but only 
the higher frequency solution with the weaker damping 
corresponds to the oscillating mode in Fig. H. 
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FIG. 4. Graphical solution of the collective mode equa- 
tion (n2|) for three different exciton densities of n — 0.05, 0.1 
and 0.2 for the parameters of Fig. pi The fine-dotted line 
corresponds to cu/n. 

The imaginary part of £ represents the damping of 
the collective mode by the two-exciton continuum. It 
peaks around w = 2uj x where the real part rises steeply 
and, thus, damps out the lower frequency solution. We 
have not included the damping by two excitons when 
one (or both) of them is dissociated into an electron-hole 
pair in our computation but this would not qualitatively 
change the damping continuum. In the three-dimensional 
model where T x ss 1 ps, these oscillations have periods 
between 50-200 fs. The curvature in the a = 0.5 case 
is exceptional, resulting in only one high-frequency solu- 
tion whose damping goes down with increasing excitation 
density. 

In conclusion, we extended the Gross-Pitaevskii equa- 
tion for the exciton polarization at low densities to in- 
clude correlation, which is an important counterpart to 
the mean-field exchange term and which provides an in- 
trinsic dephasing mechanism. A consequence is the exis- 
tence of collective oscillations riding on a slowly varying 
density. The Goldstone mode, though undetectable, in- 
dicates a Bose-Einstein condensate behavior, induced by 
optical excitation but sustained by exciton interaction. 

Finding by nonlinear optical measurement of exciton 



(12) polarization of the high-frequency mode, predicted by the 
same dispersion relation as the Goldstone mode, could 
then be argued to present evidence of condensation. 
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